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The construction of a multi-channel vector spin polarimeter requires the development of a detector 
type, which works as a spin polarizing mirror with high reflectivity and asymmetry properties to 
guarantee for a high figure of merit. Technical realizations are found by spin polarized electron 
scattering from a surface at low energies. A very promising candidate for such a detector suitable 
material consists of an oxygen passivated iron surface, as for example a Fe(001)-p(lxl)-0 surface. We 
investigate in detail the electronic structure of this adsorbate system and calculate the corresponding 
spin-polarized low-energy electron scattering. Our theoretical study is based on the fully relativistic 
SPRKKR-method in the framework of density functional theory. Furthermore, we use the local spin- 
density approximation in combination with dynamical mean field theory to determine the electronic 
structure of Fe(001)-p(lxl)-0 and demonstrate that a significant impact of correlation effects occurs 
in the calculated figure of merit. 


I. INTRODUCTION 


The first quantitative theoretical description of rel¬ 
ativistic spin-polarized low-energy electron diffraction 
(SPLEED) had been developed by Eeder (1983) and by 
Tamura and Eeder (1984) [i|,[3- The advantage of a spin- 
polarized relativistic formulation consists in the fact that 
the interplay of exchange interaction and spin-orbit cou¬ 
pling is considered on the same level of accuracy [^, Q . 
The application of this method to electron scattering 
from solid surfaces allows to support the development 
of a multi-channel vector spin polarimeter which can be 
realized using selected surfaces as two-dimensional reflec¬ 
tion mirrors. The concept of such a detection method us¬ 
ing spin-dependent electron scattering has been demon¬ 
strated, for example, for W(IOO) [5[. While in the cited 
work the detection of only one spin component has been 
realized the new scattering mirror should give the pos¬ 
sibility to detect all three spin components in a single 
step. Therefore, an optimization of suitable materials 
for the use as reflection mirrors and corresponding inves¬ 
tigations of new materials is highly desirable. The actual 
research activities focus on two classes of single crystal 
surfaces: non-magnetic surfaces from high Z-materials, 
where spin-orbit coupling acts as the underlying physi¬ 
cal mechanism or magnetic surfaces (ferromagnetic ma¬ 
terials) where both exchange interaction and spin-orbit 
coupling influence the spin-dependent electron scatter¬ 
ing. The classical system representing the spin-orbit case 
is W(IOO) Another promising candidate was found in 
the Ir(OOl) surface, which is less reactive than tungsten 
and, as a consequence, provides a longer operation time 
for the use as a spin detector Q . In contrast the scatter¬ 
ing at ferromagnetic surfaces has been investigated only 
recently [8[. Nevertheless, it was shown that using fer¬ 
romagnetic materials as spin-detectors a very high figure 


of merit (EOM) can be achieved [9[. 

In the last decades significant theoretical and experi¬ 
mental progress was made in the application of such sys¬ 
tems. One major success of these investigations was the 
determination of surface magnetic moments Q. Eor the 
exchange scattering of electrons from a sample surface 
different ferromagnetic materials have been used SET 
. An often mentioned problem of spin detectors which 
depend on spin-orbit interaction is the low EOM in the 
order of 10“^ 0 . A higher EOM was reported for ex¬ 
change scattering from an iron surface where values up to 
20 times larger have been reached [I, [TqI - [T^ . A disadvan¬ 
tage of these surfaces is their short operation time due to 
contamination. A solution to this problem is the prepara¬ 
tion of an oxygen overlayer, i.e. a surface passivation. Eor 
the coverage of 1 monolayer oxygen an ordered overlayer 
is formed resulting in p(lxl) FEED reflection patterns, 
leading therefore to a longer operation time in vacuum 
Hi- In this theoretical study the Ee(001)-p(lxl)-0 sur¬ 
face serves as a benchmark for further developments of 
our theoretical approach, as well as an suitable starting 
point for research activities on various materials which 
may applicable as reflection mirrors for spin filtering. 

The paper is organized as follows: In Sec. [TTl we de¬ 
scribe the theory of our SPLEED calculations. In Sec. 
imiwe discuss our theoretical method concerning the elec¬ 
tronic structure and the various SPLEED calculations 
and in Sec. m we summarize our results. 


II. THEORY 

A. Electronic structure 

The calculation of the electronic structure has been 
done using the Munich SPRKKR program package ET 


2 


[T^ . The implementation of the tight-binding (TB) KKR- 
method allows an effective treatment of two-dimensional 
surfaces, i.e. the self-consistent calculation of the elec¬ 
tronic structure, due to the fast convergence of the TB 
structure constants These decay expo¬ 

nentially which allows in particular the treatment of var¬ 
ious layered systems and relaxed surfaces with adsorbed 
atoms [ 2 ^. Using this method we construct a semi¬ 
infinite system with two-dimensional periodicity which 
consists of three parts: substrate (having bulk poten¬ 
tial), surface region and vacuum region (represented by 
empty spheres). The calculations were done fully rela¬ 
tivistic to treat effects coming from spin-orbit coupling 
and exchange interaction in a coherent way. To account 
for many body effects beyond the local spin density ap¬ 
proximation |2ll423j a site diagonal, nonlocal, complex 
and energy dependent self energy S determined within 
the dynamical mean field approach (DMFT) [23, [HI 
has been used. It has been shown that this method is 
straightforwardly applicable to semi-infinite lattices with 
lateral translational invariance and an arbitrary number 
of atoms per unit cell [HI, [13 • The inclusion of many- 
body effects expressed by the DMFT has been shown to 
result in significant changes for the shape of the calcu¬ 
lated spectra especially for lower kinetic energy of the 
reflected electrons [13 • An important parameter for var¬ 
ious spectroscopic calculations is the work function. We 
calculated the work function of the 2D semi-infinity sys¬ 
tem applying a summation over the Madelung potentials 
in the interaction zone as described elsewhere [HI . 


B. Theory of SPLEED 

The different asymmetries that characterize a 
SPLEED spectrum of a ferromagnetic surface are 
determined by changing the magnetization direction 
parallel to the surface either in the scattering plane 
or perpendicular to it. The scattering plane is defined 
by the wave vector of the incident and the scattered 
electrons. Beside the change of the magnetization the 
polarization of the incident electrons has to be chan^d 
resulting in four different scattered intensities 0, 3- 
These are determined by the electron polarization (a) 
as well as the direction of the magnetization (ju). It has 
been shown that mainly two different set-ups have to be 
considered concerning the orientation of magnetization 
and polarization with respect to the scattering plane 
[3- Eor our calculations we used the setup for which 
both magnetization and polarization are parallel to 
the scattering plane. If this plane is parallel to a 
mirror plane the spin-orbit asymmetry vanishes and 
the scatteriM of the electrons is only due to exchange 
interaction Q- According to symmetry considerations 
for the scattered electron intensities 


ja ^ j-a 
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holds 0- For the spin-orbit asymmetry (Agoc) defined 
by i 


Asoc — 2 ^^+ ^ 

with the definition for A_|_ and A- 
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Asoc = 0 results. The exchange asymmetry (Aex) in 
turn can be expressed by 


Aex = \{A+ + A.). (5) 

Based on the symmetry restriction in Eq. o one can 
evaluate Aex from the following simplified equation 
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As a consequence, the scattering plane is parallel to 
a mirror plane only one magnetization direction has to 
be considered when determining the exchange asymme¬ 
try. Nevertheless a useful test is to consider in addition 
the reversed magnetization to verify vanishing spin-orbit 
asymmetry. 

Another quantity to characterize different working 
points or regimes for surfaces used as scattering mirror 
is the EOM. It is defined as the product of the reflected 
intensity and the asymmetry for a specific orientation of 
the magnetization: 


FOM+(_) =/+(_). 4(_). (7) 

Here the indices indicate the magnetization direction 
of the sample. Eor the use as a spin-filter both reflectivity 
and asymmetry should have high values, leading to a high 
EOM. 


III. DISCUSSION 
A. Electronic structure calculation 

The calculation of the electronic structure has been 
started from a fully relaxed surface and interface using 
experimental structure parameters [HI- For surface sen¬ 
sitive methods it is important to include the structural 
relaxation of the topmost surface layers. Especially for 
methods using low energetic particles the changes in the 
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electronic structure resulting from the relaxations are im¬ 
portant. We set up a two-dimensional surface system 
which consists of 10 monolayers (ML) Fe, 1 ML O and 9 
ML empty spheres to simulate the vacuum. In terms of 
the TB-KKR for 2D systems we introduce left and right 
bulk regions representing the properties of the Fe sub¬ 
strate and the vacuum. For the left bulk region we used 
2 ML of Fe repeated to the left while 2 ML of the top¬ 
most empty spheres have been used for simulating the 
right bulk region, i.e. the vacuum region. Considering 
the electronic structure the interaction zone in between 
simulates the transition from Fe bulk to surface proper¬ 
ties. Our calculations were done in the atomic sphere 
approximation (ASA) using the parameterization for the 
exch ang e-correlation functional of Vosko, Wilk and Nu- 
sair j^. We used a lattice constant of 2.86 A according 
to the unit cell of bee Fe. A fully relativistic ab initio 
calculation of the potentials was applied to account for 
spin-orbit and exchange effects on one footing. We also 
included many body effects in our calculations consider¬ 
ing the sensitivity for spectroscopies based on low ener¬ 
getic particles. For solving the many-body effective im¬ 
purity problem a DMFT solver has been used [3l|, [32| . In 
our calculations we utilized the spin-polarized T-matrix 
approximation solver (TMA) [^. It has been shown 
that the use of a TMA solver is justified because of the 
less pronounced correlation effects in transition metals 
[33 • The parameters which have to be supplied for a 
DMFT calculation is the intra-atomic Hund exchange in¬ 
teraction (J) and the screened Coulomb interaction (U). 
Corresponding to previous extensive studies we set the 
values to J = 0.9 eV and 1/ = 2.3 eV j^ . 

The results of the electronic structure calculations are 
shown in Fig. [T]in terms of the density of states (DOS) of 
the first three Fe layers and the O layer. Beside the result 
of LDA calculation the DOS of a DMFT calculation is 
shown as well. 

The DOS shows reasonable agreement with previous 
electronic structure calculations for the system Fe(OOl)- 
p(lxl)-0 [35|. It is visible that in the energies regime 
—6 eV to —2 eV a hybridization of the valence states 
between O and the Fe layer occurs. The LSDA+DMFT 
based DOS calculations show a broadening especially for 
the topmost Fe layers resulting in spectral changes for 
surface sensitive spectroscopic methods. This is caused 
by the finite value of the self energy in the specific energy 
range and is also visible in the calculated band struc¬ 
ture (see Fig. [8] and Fig. [9]). The changes in both 
spin channels result in a lower spin magnetic moment 
{i^spin) for the DMFT calculations. The main features 
agree using LSDA or LSDA+DMFT indicating that the 
main properties of the electronic and magnetic structure 
of the Fe(001)-p(lxl)-0 surface can be described using 
both schemes. Nevertheless, for the calculation of very 
low-energy electron diffraction it is important to include 
LSDA+DMFT (see below). In Fig. [2] the spin- and or¬ 
bital magnetic moments for bee Fe bulk and the three 
topmost atomic layers of Fe(001)-p(lxl)-0 are shown. 
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FIG. 1. DOS of the first three atomic layers of the inves¬ 
tigated half-infinity surface system Fe(001)-p(lx)-O. In a) 
the DOS for a standard LSDA calculation is shown. In b) 
the calculation includes many-body effects accounting via 
LSDA-hDMFT. The DOS of O and the first Fe layer are drawn 
together in the uppermost panel to show the hybridization of 
O and Fe in the valence band. The inset shows a sketch of 
the surface system and the numbering of the atomic layers. 


The Fe atoms of the topmost layer have a higher magnetic 
moment when compared to the bulk value. This finding 
is known for magnetic atoms on surfaces and is related 
to the band narrowing of the d-states [36|, [13 • The trend 
for the decrease of the spin magnetic moments going to 
deeper Fe layers is similar for LSDA and LSDA-f DMFT 
calculations. This behavior is reflected by the compa¬ 
rable main features in the DOS (see Fig. [1]) comparing 
LSDA and LSDA-fDMFT. 
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FIG. 2. Magnetic moments for Fe bulk and the first three 
atomic layers of the semi-infinite system Fe(001)-p(lxl)-0 
with and without the inclusion of many-body effects. Beside 
the spin magnetic moments the orbital moments are shown. 
For O the orbital moment has zero value for LDA and DMFT. 


Due to the hybridization of O and Fe in the valence 
band regime a magnetic moment for O is induced. The 
increased magnetic moment of Fe for the surface layer in 
combination with the induced magnetic moment for O 
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results in a larger exchange interaction at the passivated 
Fe(OOl) surface in contrast to a clean Fe(OOl) surface. 
Comparing the spin-orbit induced orbital magnetic mo¬ 
ments the differences between LSDA and LSDA-f-DMFT 
calculations are more pronounced than for the spin mo¬ 
ments. Using the LSDA an enhancement of the orbital 
moment for the two outermost atomic layers occurs as 
one would expect. Beside the decrease of the orbital mo¬ 
ment going to deeper Fe layers is stronger for the LSDA 
calculation. Nevertheless the dominating part is the spin 
magnetic moment which characterises the exchange scat¬ 
tering of the polarized electrons. In summary the passi¬ 
vation of the Fe surface results in a significant change of 
the magnetic properties compared to a non-passivated Fe 
surface. In our calculations for the first Fe layer of a non- 
passivated Fe surface a spin and orbital magnetic moment 
rris = 2.S1 fiB and rrio = 0.11 fiB results. Whereas for a O 
passivated surface using the LSDA for the first Fe layer 
rris = 3.25/iB and rrio = 0.01 jiB have been calculated. 
The enlarged magnetic moments results in an increased 
exchange scattering as been shown in previous experi¬ 
ments [1^. Also a very high Sherman function was re¬ 
ported which is linked directly to the magnetic properties 
of the surface [9[. It should be mentioned that our cal¬ 
culated magnetic moments especially for the two surface 
layers show good agreement with data in the literature 

[H. 


B. SPLEED calculation 

For the SPLEED calculations, the scattering plane was 
aligned along the [100] direction whereas the surface mag¬ 
netization as well as the polarization of the electron was 
aligned along the [±100] direction. All calculations were 
done for the specularly reflected beam, i.e. the (0,0) 
beam using the the surface potential barrier of Rundgren- 
Malmstrom [^. Eor the O passivated Ee(OOl) surface we 
calculated a work function of 7.07 eV. In comparison to 
a clean Ee(OOl) surface a value of 5 eV was calculated, 
i.e. an increase of the work function by passivation was 
found. It should be mentioned that the work function of 
Ee(OOl) shows reasonable agreement with experimental 
and theoretical values in the literature [39[. We calcu¬ 
lated polar angle {0) - energy maps that are shown in 
Fig. [3]- Eig. [5] for the reflectivity, the exchange asym¬ 
metry and the EOM, respectively. The polar angle was 
varied in the range from 27-75° whereas the energy range 
was set to 1.3 to 17 eV according to the possible working 
areas as scattering mirror. 

In the right panel of Eig. [3] the effective reflectivity 
of Ee(001)-p(lxl)-0 is shown. As can be seen we get a 
huge reflectivity especially for relatively low kinetic en¬ 
ergies over the full range of polar angles. In particular 
at a kinetic energy of 6 eV and a polar angle of 30° a 
maximum of the reflectivity occurs. Also visible in Eig. 
[3] is the emergence threshold starting around 4 eV and 
a polar angle of 75° and ending at 8 eV and 27° which 


marks the occurrence of a new beam. It divides the map 
into mainly two parts of higher and lower reflectivity. 
This results from the fact that for kinetic energies above 
the emergence threshold the additional scattering chan¬ 
nel lowers the intensity of the specular beam as shown in 
recent experiments [^. The left panel in Eig. [3] shows 
the effective reflectivity for a clean Ee(OOl) surface. Eor 
the passivated Ee surface higher values for kinetic en¬ 
ergies greater than the emergence threshold have been 
obtained. This is due to the higher exchange scattering 
for the oxygen passivated Ee surface coming out of the 
higher magnetic moments at the topmost atomic layers. 
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FIG. 3. Left: ©-energy-map of the reflectivity for Fe(OOl). 
Right: The same for Fe(001)-p(lxl)-O. 


In Fig. 0] the exchange asymmetry for Fe(OOl) and 
Fe(001)-p(lxl)-O, respectively, is shown. The plots in¬ 
clude different areas according to the preferred reflected 
spin-orientation for a defined orientation of the magneti¬ 
zation. 
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FIG. 4. Left: Exchange asymmetry (A+) for Fe(OOl). Right: 
The same for Fe(001)-p(lxl)-O. 

The positive values of the exchange asymmetry corre¬ 
spond to a parallel alignment of the electron spin and the 
sample magnetization. Eor Ee(001)-p(lxl)-0 in the en¬ 
ergy region below the emergence threshold the scattering 
of parallel aligned electron spin and magnetization is pre¬ 
ferred except for a small area located at the emergence 
threshold. Eor polar angles larger then 50° a crossing 
of the emergence threshold results in a change of the 
scattering behaviour. In this case the reflected spin di¬ 
rection can be rotated by changing the kinetic energy of 
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the diffracted electrons. For kinetic energies higher than 
6 eV the scattering of electrons with antiparallel spin 
alignment is preferred. Beside of a change of the polar¬ 
ization of the electron a change of the magnitude occurs. 
The values of the calculated exchange asymmetry fit well 
to the experimental data [40| . 

Comparing the results for the asymmetry of Fe(OOl) 
and Fe(001)-p(lxl)-0 one notes that for the complete 
range of energy and polar angles the asymmetry changes. 
This is due to the different magnetic properties of the 
Fe(OOl) and the Fe(001)-p(lxl)-0 surfaces. At kinetic 
energies below the emergence threshold the asymmetry 
for Fe(001)-p(lxl)-0 shows broader areas with one spe¬ 
cific asymmetry direction. This is in line with to other 
investigations made for these systems 0- 

In Fig. [5] the FOM is shown which is the most impor¬ 
tant observable for characterizing a material to be used as 
possible spin-filter. On the left side the FOM for Fe(OOl) 
is shown whereas on the right side for Fe(001)-p(lxl)- 
O. For the oxygen passivated surface a broad range of a 
high FOM for kinetic energies lower than the emergence 
threshold occurs. This comes from the asymmetry A+ 
(see Fig. 0]) according to Eq. ([7]). It is an advantage 
for the application as spin filter using low kinetic en¬ 
ergy electron diffraction for the determination of surface 
properties. For Fe(001)-p(lxl)-0 the highest values are 
reached for kinetic energies lower than 4 eV and a po¬ 
lar angle larger than 50°. Due to the fact that working 
points for spin-filters are suitable between 40° and 60° 
the oxygen passivated surface is a promising candidate 
for such an application. 
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FIG. 5. Left: ©-energy-map of the exchange FOM for 
Fe(OOl). Right: The same for Fe(001)-p(lxl)-O. 

In Figs. [6] and [7| the experimental results for the reflec¬ 
tivity, the FOM and the exchange asymmetry are shown 
for comparison jdoj . 

The emergence threshold is well reproduced by the the¬ 
oretical results, i.e. its correspondence for varying polar 
angle and kinetic energy. Hence the geometric configura¬ 
tion described by our calculations match the setup in the 
experiment. Also the inner potential calculated out of the 
work function and the Fermi energy is confirmed. This 
is ensured by the fact that a difference to the experiment 
would result in a energy shift of the emergence threshold. 
In Fig. [7] the change of the exchange asymmetry is shown 
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FIG. 6. Experimental results for the SPLEED measurements 
on Fe(100)-p(lxl)-0 taken from Ref. 0- (Reproduced by 
permission.) Left: ©-energy-map of the reflectivity. Right: 
©-energy-map of the exchange FOM. 


by reversing the magnetization of the Fe(100)-p(lxl)-0 
surface. The same behavior as in the theoretical results 
are visible coming out of vanishing spin-orbit asymmetry. 
Based on that changing the magnetization direction ex¬ 
actly inverses the scattered polarization of the electrons. 
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FIG. 7. Experimental results for the SPLEED measurements 
on Fe(100)-p(lxl)-0 taken from Ref. (Reproduced by 

permission.) Left: ©-energy map of the exchange asymme¬ 
try. Right: ©-energy map of the exchange asymmetry for the 
reversed magnetization. 


C. Many-body effects in SPLEED calculations 

The importance of the inclusion of many body effects 
for spectroscopic calculations has been shown in differ¬ 
ent works [26[ 113, l4ll - [43| . Because of the low kinetic 
energy of the incident electrons many-body effects might 
become important for the spectroscopic calculations due 
to changes of relevant bands (surface states, bands at the 
Eermi energy, unoccupied states) resulting from a change 
of the underlying electronic structure calculations. In 
Fig. [8] the band structure of Ee(001)-p(lxl)-0 with and 
without inclusion of the DMET is shown. As one notes, 
the bands around 7 eV above the Eermi level are smeared 
out by many body interactions. This energy range is im¬ 
portant for SPLEED giving characteristic spectral fea¬ 
tures. Changes in the band dispersion results in differ- 
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ent magnetic properties altering the exchange interac¬ 
tion at the sample surface. This effect can be connected 
to the self energy shown in Fig. [9l Comparing to the 
calculated work function (7.07 eV) the self energy has 
a nonzero value affecting the valence bands relevant for 
the exchange scattering process. Although our calcula¬ 
tions show that the impact on the effective reflectivity 
and the exchange asymmetry is negligible, the changes 
in the FOM are significant. 
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comparison to Fig. [5]this results in a shift of the maximal 
value of the FOM to a polar angle around 50°. Accord¬ 
ing to the changes of the band structure which affects 
essentially the band near the Fermi level, changes in the 
FOM are seen mainly for low kinetic energies. This is 
important comparing measurements and calculations of 
spectra for low energy electron diffraction. 
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FIG. 10. FOM without (left) and with (right) consideration 
of many body effects for Fe(001)-p(lxl)-O. In both cases the 
FIG. 8. Left: Bloch spectral function without consideration projection of the electron spin has been included, 
of many body effects for Fe(001)-p(lxl)-O. Right: Results in¬ 
cluding many body effects via the DMFT for Fe(001)-p(lxl)- 
O. 

IV. SUMMARY 
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FIG. 9. Real and imaginary part of the self energy for both 
spin channels. 


In Fig. [To] we present the FOM for Fe(OOl)- 
p(lxl)-0 resulting from LSDA+DMFT-based SPLEED- 
calculations. Additionally we considered the change in 
the projection of the polarization of the electron con¬ 
cerning the surface when changing the polar angle. In 


We have shown that the calculations done using our 
ab initio method regarding the SPLEED-spectra for 
Ee(001)-p(lxl)-0 are in satisfying agreement with recent 
experimental results [40|. Therefore our description of 
the systems electronic properties seem to be confirmed. 
The system exhibits a large EOM and various suitable 
areas for the application as spin-polarizing mirror. We 
have shown that a projection of the polarization of the 
electron has a huge impact on the exchange scattering, 
especially for the calculation of the EOM. Eurthermore 
the inclusion of many body effects has been considered 
for the first time in SPLEED calculations showing that 
this results in spectral changes important for the regime 
of very low energy electron scattering. 
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